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ABSTRACT 

Spectral observations of molecular line profiles reveal the so-called 'blue profiles' for 
double-peaked molecular lines with stronger blue and weaker red peaks as notable 
features for star-forming cloud core collapses under the self-gravity. In contrast, 
25 — 30 per cent of observed molecular spectral line profiles in star-forming clouds 
or cores also show the so-called double-peaked 'red profiles' with red peaks stronger 
than blue peaks. Gao & Lou (2010) show that these unexplained 'red profiles' can 
be signatures of global dynamics for envelope expansion with core collapse (EECC) 
within star-forming molecular clouds or cores. We demonstrate here that spatially- 
resolved 'red profiles' of HCO+ (J= 1 — 0) and CS (J= 2 — 1) molecular transitions 
from the low-mass star-forming cloud core FeSt 1-457 together with its radial profile 
of column density inferred from dust extinction observations appear to reveal a self- 
similar hydrodynamic shock phase for global EECC. Observed spectral profiles of 
C^^O (J= 1 — 0) are also fitted by the same EECC model. For further observational 
tests, the spatially-resolved profiles of molecular transitions HCO+ (J=: 3 — 2) and 
CS (J= 3 — 2) as well as the radial profiles of (sub)millimetre continuum emissions at 
three wavelengths of 1.2mm, 0.85mm and 0.45mm from FeSt 1-457 are also predicted. 

Key words: ISM: clouds — ISM: individual (FeSt 1-457) — line: profiles — radiative 
transfer — stars: formation — stars: winds, outfiows 



1 INTRODUCTION 

Gravitational collapses of molecular clouds and cloud cores 



lead to early formation of pr o tostars therein (e.g. Shul 



Shu, Adams fc Lizanol 1 19871 : iMcLaughlin fc Pudrit? 



1977 



1997 



McKee fc Tanll2002l ). For starless and star-forming molec- 



ular cloud cores in early phases, the central dense mass 
blob might not involve discs or outflows especially for those 
grossly spherical molecular clouds. Recent theoretical de- 
velopment of self-similar hydrodynamic shock models in- 
volving the self-gravity reveals, among others, a novel plau- 
sible dynamic scenari o of global envelope expansion with 
core colla pse (EECC) (Lou fc Shcn"2004|; lBian fc LoulbOOSi : 
ILou fc Gao 2006 : Wang fc Lou 2008) with an expanding in- 
fall radius Rint during early dynamic evolution stages of 
forming protostellar cores. Together with other pertinent 
yet independent diagnostics and constraints, observations of 
various molecular emission spectral line proflles, especially 
those with high enough spatial and frequency resolutions 
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and/or different line transitions of the same molecules, make 
it possible to identify the gross hydrodynamic and thermo- 
dynamic structures of star-forming molecular clouds with 
grossly spherical morphologies. 

Major diagnostics to probe molecular cloud structures 
include spatially-resolved s pectral profiles for various molec- 
ular line transitions (e.g. IZhou et al.l Il993l : iTafalla et al.l 
|2006| . Fu, Gao & Lou 2010), radial profiles of column densi- 
ties inferred from dust extinctions of backgr ound stars be- 
hind a molecular cloud (e.g. Alves et al. 2001: iKandori et al.l 
I2OO5I), radial profiles of ( s ub)millimetre continu um emissions 
fe.g. IShirlev et aLlbOOol : iMotte fc Andrei 120011 ) etc. Any of 
such separate observational diagnostics may be fitted with 
empirical and parametrized models on an individual basis 
to various extents as can be found in the literature. In con- 
trast and in order to remove potential model degeneracies, 
our strategy is to constrain the underlying cloud model us- 
ing the hydrodynamics and by all available observational 
data simultaneously. We also make model predictions for 
specific observational tests. We mainly focus on those molec- 
ular clouds which appear grossly spherical as our spheri- 
cally symmetric hydrodynamics represents a first-order ap- 
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proximation in that direction. For illustration, we choose 
the molecular cloud core FeSt 1-457 for specific observa- 
tional comparisons. We also summarize the main results of 
our model comparisons with several observations for another 
star-forming cloud core L1517B as supporting evidence. 

Regarding the molecular line diagnostics, there are two 
major asymmetric double-peak molecular spectral line pro- 
files oft-observed in star-forming clouds, viz., the so-called 
'blue profiles' with stronger blue-shifted peaks than red- 
shifted peaks and the so-called 'red profiles' with stronger 
red-shifted peaks than blue-shifted peaks. The level of 
such spectral line profile asymmetries can change in a 
gradual manner for spatially-resolved observations. While 
the observed 'blue profiles' of molecular spectral lines are 
modelled and generally interpreted as ch aracteristic signa- 
tures of r adial infalls in molecular clouds (|Zhou et al.lll993 : 
iGao. Lou fc Wu 2009), the somewhat ignored yet wide pres- 
ence of unexplained 'red profiles' (e.g. iMardones et al.iri997l : 
lEvansll2003l : I Velusamv et al.|[2008l ) remains a conceptual gap 
in our understanding of star-forming molecular cloud dy- 
namics. Moreover, sometimes 'blue profiles' and 'red pro- 
files' coexist i n molecular cloud cores for the same line 
transition (e.g. iTafalla et al]l2000l). and/or for different line 
transitions (e.g. iTafalla et all 20061 ). For the possible inter- 
pretation of 'red profiles' of molecular spectral lines alone, 
one may thin k of several mechanisms, e.g. rotations (e.g. 
iDi Francesco et al. 2001i '). b ipolar outfiows inclined c lose to 
the plane of the sky (e.g. Fiege &: Henrikse"nl llOgd ). tur- 
bulence (e.g. | Ossen kopfll200i r and thermal pulsations (e.g. 
iLada et al.ll2003l l. 

For FeSt 1-457 (see Fig. 1), disc rotation is unlikely, as 
this cloud core is most likely evolving in an early phase. 
A global cloud core rotation of FeSt 1-457 may be ex- 
cluded by spatially-resolved spectral line observations; i.e., 
the currently data do not show such evidence. At present, 
these mechanisms just mentioned have not been shown to 
also account for cloud density profiles and radial profiles of 
(sub)millimetre continuum emissions simultaneously within 
a consistent model framework. Alternatively, it was sug- 
gested that 'red profiles' in certain star- for ming sources 
may be caused by clouds invo lving EECC (jLou fc ShenI 
l2004l : iThompson fc White! liooi ) . It is thus crucial to re- 
produce the observed 'red profiles' using our general poly- 
tropic EECC_slwck_jnodels via radiative transfer calcula- 
tions iGao fc Loull2010l ). Meanwhile, we require the model 
density profile to be consistent with the observational infer- 
ence of dust extinctions and further predict other spectral 
line profiles as well as radial profiles of (sub)millimetre con- 
tinuum emissions as future observational tests. 

Four main contributions of our model analysis in this 
paper are as follows. First, for the grossly spherical molecu- 
lar cloud core FoSt 1-457, we construct a self-similar gen eral 
polyt ropic liydrodynami c EECC shock mod el (Lou fc S henI 
120041 : (Lou fc Gadl2006l : I Wang fc Loi]l2008l L and radiative 
transfer calculations are performed to reproduce the ob- 
served 'red profiles' and other spectral profiles of molecular 
transition lines in the cloud core. Secondly, we advance the 
conceptual framework of constraining the molecular cloud 
dynamics using several available diagnostics simultaneously 
in order to solve this challenging issue in molecular astron- 
omy. For FeSt 1-457, the column density profile of our EECC 
shock dynamics is also sensibly fitted with the observa- 



tionally inferre d column density pr ofile by the dust extinc- 
tion data (e.g. iKandori et al]|2005h . Thirdly, while no data 
for (sub)millimetre continua are available for FeSt 1-457 at 
present, for scenario consistency and further observational 
tests, we predict radio continuum emissions at three wave- 
lengths of 1.2 mm, 0.85 mm and 0.45 mm. For the same 
reason, the spatially-resolved spectral profiles of molecular 
line transitions HCO+ (J= 3-2) and CS (J= 3-2) are 
also predicted for observational tests. Finally, we suggest 
(|Gao fc Loull2010l ) that such EECC shock phase may evolve 
from damped nonlinear radial pulsations of molecular cloud 
core within one or two cycles. Corresponding to different 
phases of such damped large-scale radial pulsations, differ- 
ent dynamic phases may give rise to distinct diagnostics. 



2 STAR-FORMING MOLECULAR 
CLOUD CORE FEST 1-457 

As a chosen candidate, the molecular cloud core FeSt 1-457 
appears to be a grossly spherical dark globule in the Pipe 
Nebula at a distance of ~ 130 pc away (e.g. Lombardi et 
al. 2006). Without detecting a star, this starless cloud core 
in the early phase of protostar formation has a diameter of 
~ 35, 000 au with an angular diameter ~ 4.5' derived from 



the o ptical images (e.g. l Aguti et al..,2007l') . The spectral data 



3pt 

from lAguti et al.l (120071 ') for spatially-resolved spectral line 



profiles of three molecular transitions HCO (J= 1 — 0), 
CS (J= 2-1) and C^**0 (J= 1 - 0) (Figs. [1] and g]) were 
acquired in 2003 by the 30-metre millimetre-wave telescope 
of Institut de Radioastronomie Millimetrique (IRAM) on 
Pico Veleta near Granada in Spain. The spatial resolution 
of this spectral data set is ~ 24" for the beam width in the 
frequency band 90 — 110 GHz with a mean beam efficiency of 
~ 0.75 at uncertainties < 4 per cent (e.g. Greve et al. 1998). 
The dust extinction data of FeSt 1-457 for the inferred radial 
profile of column density with a ^ 30" beam resolution are 
also available together with error bars (e.g. iKandori et al.l 
[2005i . and see the upper panel of our Fig. (HJ. 

Global expansions in the outer layer or envelope of this 
globule are hinted by asymmetric molecular spectral line 
profiles and cloud oscillatio ns were proposed to be the ori- 
gin of such expansions (e.g. lAguti et al.ll2007h . In our advo- 
cated scenario, such large-scale nonlinear cloud radial pulsa- 
tions cannot persist forever but may initiate core collapses 
and envelope expansions within one or two cycle s and hence 
may l ead to a self-simil ar EECC shock phase (|Lou fc Gaol 
'2006?; iGao fc Loull2010l l. We demonstrate that the global 
EECC shock phase seem to persist in FeSt 1-457 by molec- 
ular spectral line profile fits using the spherica lly symmetric 
;eneral polytropic self-similar dynamic model ijWang fc Lou! 
2008 ). Meanwhile, w e also use the dust extinction data 



Kandori et al.ll2005l l of FeSt 1-457 for the radial profile of 
column density with a ~ 30" beam resolution to constrain 
our EECC shock hydrodynamics (Fig. |6]). Finally for future 
observational tests, we predict the spatially-resolved spec- 
tral profiles of molecular line transitions HCO"'' (J= 3 — 2) 
and CS (J= 3—2) as well as three radial profiles of millimetre 
continua at wavelengths 1.2mm, 0.85mm and 0.45mm with 
estimated error ranges using the density and temperature 
profiles of the same underlying EECC shock dynamics (Fig. 
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Figure 1. The observed molecular line transition HCO"*" (J= 1 — 0) in a spatially-resolved grid- map of star- forming molecular cloud 
core FeSt 1-457 (solid line profiles) superposed with our best-fit EECC shock dynamic model results (dotted line profiles). The receding 
velocity of the entire cloud core FeSt 1-457 along the line of sight is estimated to be ~ 5.75 km s~^. Each grid has the local receding 
velocity along the line of sight as the absciss a and antenna temperature as the ordinate (see the very upper right grid as a sample). The 
observational data map of lAguti et al.l l|2007l ) is centered at right ascension « j2000 = 17''35™47.5= and declination 5 

.J2000 = — 25°33'2.0", 

with a grid spacing of ~ 24" (corresponding to ~ 3.0 X 10^ au at a distance of ~ 130 pc for the Pipe Nebula) and a maximum observed 
area of 168" x 168" in the plane of the sky . Molecular spectral line profile fittings are based on our general polytropic EECC shock 
model llWane &: Loull2008l : iGao fc Loull2010l) . by adopting a step function molecular abundance profile (see details in Table 1) and by 
performing the RATRAN computation. Our model spectral line profiles are convolved with the beam width of 24" and multiplied by 
an IRAM 30m telescope mean beam efficiency of ~ 0.75. 



[5]and|6]), as have been already det ected in other proto star- 
forming molecular cloud cores (e.g. IShirley et al.ll2002D . 



3 EECC SHOCK HYDRODYNAMICS AND 
FORMATION OF RED PROFILES FOR 
MOLECULAR SPECTRAL LINES 

3.1 Self-similar Hydrodynamic Models 
and EECC Shock Solutions 

EECC self-similar hydrodynamic shock solutions have 
been investigated with different equa ti ons of state (EoS 
in recent years JLou fc ShenI |2004 ILou fc Gaol 1200' 
IWang fc Loull2008h . We here adapt the magnetohydrody- 
namic (MHD) mode l form ulation with general polytropic 
EoS in IWang fc Loul (120081 ) but without the magnetic field 
(i.e. setting the magnetic field parameter /i = in that pa- 
per). 

Our model formulation begins with the ideal hydrody- 
namic partial differential equations (PDEs) for a spherically 
symmetric molecular cloud in spherical polar coordinates 
(r, 9, <j)), namely 



21 + 1. 

dt r2 



(1) 



du du 
dt dr 



Idp 
p dr 



GM 



dM ^ dM _ ^ 
dt dr 



dM , 2 
— = 4^r p , 



(2) 
(3) 
(4) 



where, as functions of both radius r and time t, dependent 
variables u, p, M, and p are the radial flow velocity, mass 
density, enclosed mass within radius r at time t and ther- 
mal gas pressure, respectively; G is the gravitational con- 
stant and 7 is the polytropic index. These nonlinear hy- 
drodynamic PDEs describe the conservations of mass (eq 
m), radial momentum (eq 0) and specific entropy along 
streamlines (eq Q). 

Self-similar solutions to nonlinear PDEs (H} — Q by 
the following self-similar transformation can be constructed, 
namely 



u = k^'h 



a(x) 



(5) 



(6) 
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M 



(3n - 2)G 



kt^''-^a{x)''m{x)'' 



(7) 



where x is the self-similar independent variable combining 
r and t in a special manner, and v{x), a{x) and m{x) are 
functions of x only and represent the dimensionless reduced 
radial flow velocity, mass density and enclosed mass, respec- 
tively; and n and 7 in these self-similar transformations are 
two additional scaling indices allowed, with the algebraic re- 
lation q — 2(n-|-7 — 2)/(3n — 2). Here k is a parameter related 
to the sound speed within the molecular cloud. Physically, 
a positive enclosed mass M requires the inequality n > 2/3 
in the similarity transformation above. For n -|- 7 = 2 and 
thus g = as in this paper unless otherwise stat ed, we reduce 
to the conv entional polytropic hydrodynamics (|Wang fc LorJ 
l2007l .[ 20081 ) as an important subcase. Other important phys- 
ical variables, namely the gas thermal temperature (invoking 
the ideal gas law) and the central mass accretion rate in the 
molecular cloud centre are 

E^kt^(--^^a{xr-'m{xy , (8) 

(9) 



kBp/ifimn) kB 
Mo = fc^'''^t^<""''mo/G 



respectively. Here mo = m(0) is the reduced central enclosed 
mass, fca is Boltzmann's constant, mn is the atomic hydro- 
gen mass, and /i is the mean molecular weight setting to be 
2 in our model analysis for star-forming gas cloud consisting 
mainly of molecular hydrogen H2. 

By adopting self-similar transformation equations 
((SI — (O, nonhnear hydrodynamic PDEs ([l]) — Q can be 
readily cast into two coupled nonlinear ordinary differential 
equations (ODEs) which can be solved by a c ombination of 
analy tical and numerical techniques (see, e.g. IWang fc LorJ 
l2008t ). These two coupled nonlinear ODEs are 



2-n + 



3n - 2 



q \a 



2 x ''{nx — v)'' — {nx — v) 



da „(x — v) 

X— = 2^ -a 

ax x 



{nx — v) + qa^ "'^ 2 x^''{nx — v 



, , > (nx — v) 



+qa 



x^'' ''(nx — vy ^(3nx — 2v) 



(10) 



3n — 2 \ 2<7/ NO / ^2 

2 — n ~\ q]ct 2 X (nx ~ v) — (nx — V 



dv Ax-v)( 3ng \ 2q, .q 

X — — 2- 2 — n H a ^ x ^{nx ~ v) 

dx X \ 2 ' 



'{nx — v) 



. ^ , (nx — v) 



+qa 



l-n+^ 29-1 



X {nx — vY (Sns — 2v) 



(11) 



In the above mathematical derivations, we also make use of 
the following algebraic relation 

m = ax^{nx — v) (12) 

from the mass conservation ([1]) and ((3]). Once a{x) and v{x) 
axe known, we then readily obtain the reduced enclosed mass 



m,{x). Two sets of analytical asymptotic similarity solutions 
for X -> +CX} and a; — >■ 0^ respectively serve as the boundary 
conditions in solving two coupled nonlinear ODEs (jlOp and 
(lll|) . Among other possible solutions, we have for the x — > 
+00 regime 

a — Ax~^^", 



nA 



+ 2(2 - 71)71"-^^-"+^"'"'^ 



(n-2)/n 



(3n - 2) 

and for the x 0^ regime, the central free-fall solution 



2mo 



(3n - 2)x 



1/2 



(3n — 2)mo 
2^3 



1/2 



(14) 



Here, the mass parameter A and the velocity parameter B, 
describing the asymptotic boundary condition as x ^ +00, 
are two constants of integration and are important parame- 
ters to characterize the global self-similar dynamic solution. 
Meanwhile, the reduced enclosed mass mo = m(0) describes 
the central point mass embedded deep in the cloud core. 
Then our global self-similar dynamic solutions can be con- 
structed by straightforward numerical integrations using the 
standard fourth-order Runge-Kutta scheme and by specify- 
ing boundary conditions p3[) and p4|) . In the absence of 
a shock, A and B are specified and mo is determined ac- 
cordingly; i.e. only two are independent among the three 
parameters. For solutions with shocks, the shock location 
is another parameter and discontinuity across a shock front 
is introduced; global solutions can be still constructed by 
simply requiring the conservation laws of mass, radial mo- 
mentum and energy across the upstream and downstream 
sides of the shock front (see Wang fc Loul (2008) for more 
details regarding the hydrodynamic formulation). 



3.2 The EECC hydrodynamic Shock Solution 
Adopted for the Cloud Core FeSt 1-457 

In our model analysis, various EECC shock dynamic solu- 
tions are extensively explored in order to plausibly represent 
the dynamic state of FeSt 1-457 by adopting model fitting 
procedures detailed in Section 4. Finally after an extensive 
search, a sample EECC shock dynamic solution with poly- 
tropic index 7 = 1.2, scaling parameter n = 0.8 and three 
coefficients of asymptotic conditions A = 4.65, B — 2.05 and 
mo — 0.038 is choserQto fit the molecular spectral line pro- 
file data and the column density radial profile inferred from 
dust extinction obser vations. Two relev ant physical scalings 
(see subsection 2.3 of ICao fc LoullioiOl ) of the length scale 
and number density scale are specifically chosen below as 



{AnGfimut'^y^ = 2.7 x 10"* cm"^ 



(15) 
(16) 



with the mean molecular weight jj, — 2. 

Our EECC hydrodynamic shock solution offers the 
global radial structural profiles of radial infall velocity 



^ The outgoing shock radius and thus shock travelling 
current epoch t are then determined. 



speed at 
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Figure 2. Radial profiles of infall speed — t), H2 number den- 
sity A''(r, i) and thermal temperature T{r, t) for the polytropic 
EECC shock dynamic phase of a model cloud for FeSt 1-457 with 
an evolutionary timescale t ~ 2.8 X 10^ yr. This EECC shock dy- 
namic model has an expanding infall radius -Rinf = 4.5 X 10^ au 
(top panel), a shock radius Rsh = 8.16 X 10^ au (all three pan- 
els) and an outer cloud boundary of -Rout = 1.80 X 10* au. 
The infall radius Ri^f and shock radius iigh expand at speeds of 
~ 0.015 km s~^ and ~ 0.27 km s~^, respectively. The top panel 
shows that the cloud core collapses inside Rinf while the enve- 
lope outside Rinf expands outwards with typical gas flow speeds 
< 0.1 km s~^. The middle and bottom panels show a number 
density increase by a factor ~ 10^ and a temperature increase by 
a factor > 3 times towards the core centre where a point mass 
blob of Mo = O.O2O2M0 resides. 




Figure 3. A schematic illustration of the envelope expansion with 
core collapse (EECC) scenario for generating red profiles of asym- 
metric molecular line transitions. The cloud consists of an outer 
expanding envelope and an inner collapsing core, separated by an 
outgoing spherical surface at r = r^^i^n in a self-similar radial ex- 
pansion. The loci (i.e. the two connected partial ovals inside the 
infall radius rinf^n) and two straight arrow lines pointing right of 
particularly chosen constant projected line-of-sight velocity com- 
ponents are shown as examples. Line of Sight 1 passes through 
the core of central infall as well as the outer expanding enve- 
lope, that intercepts the ovals and straight arrow lines at three 
points Bi, B2 and B3 for the same approaching projected velocity 
components and at three points Ri, R2 and R3 for the same re- 
ceding projected velocity components. Sufficiently far away from 
the cloud centre. Line of Sight 2 passes completely through the 
outer expanding envelope. 



u{r, t), number density N{r, t) and thermal gas temper- 
ature r(r, t) (see Fig. [2|; the cloud temperature T{r, t) 
increases towards the centre and an infall radius i?inf = 
4.5 X 10^ au expanding at a speed of ~ 0.015 km s~^ sepa- 
rates the inner core collapse domain {u < 0) and the outer 
envelope expansion domain (u > 0). An outgoing spheri- 
cal shock is at radius i?sh ~ 8.16 x 10^ au with a shock 
travelling speed of ~ 0.27 km s~^ and the outer cloud core 
radius is at iZout = 1.80 x 10'' au. Following the chosen scales 
by expressions (fT^ and (fTU)) , the sound parameter k values 
across a shock front are different; the upstream (pre-shock) 
side has k^^"^ = k]/^ = 48.7 km s"*^ * and the downstream 
(post-shock) side has fc^/^ = fc^^ = 49.1 km 8"° ^ Cor- 
respondingly (also as a result of different densities in the 
upstream and downstream sides, see eq. (8) in Lou & Gao 
2006), the sound speed c across the outgoing shock front is 
different: Cu — 0.133 km s"'^ and Cd = 0.143 km s^^. This 
EECC shock dynamic model describes a molecular cloud 
core in the early phase of star formation (a timescale of 
t ~ 2.8 X 10^ yr) with a central mass accretion rate of 
Mo ~ 0.69 X 10~^ Mq yr~^ and a central point mass blob 
of Ma = 0.0202 Mq. 

For the identification of a sensible self-similar EECC 
shock model for FeSt 1-457, we have explored different in- 
dependent parameters within the model framework, namely 
the polytropic index 7, the mass and velocity parameters 
A and B for the asymptotic boundary condition at large x, 
and a self-similar radius xi describing the shock location and 
thus the shock speed. Then the other two indices, n and 7 



are naturally fixed under the condition of conventional poly- 
tropic case with n-f 7 = 2 and q — 0; and the inner reduced 
point mass mo is also determined by the global trend of the 
specific EECC shock solution as A, B and shock location 
x-i are chosen. There are additional independent parameters 
in the scaling process of the self-similar model, namely the 
length scale (|15|) and the number density scale (|16|) . The 
current set of independent parameters for the EECC shock 
dynamic model, without degeneracy as we know, sensibly 
describes physical properties within molecular cloud core 
FeSt 1-457 when we compare the hydrodynamic shock cloud 
model with currently available observations as described in 
Section 4 presently. We further venture to make pertinent 
predictions based on our theoretical model analysis for fu- 
ture observational tests. 



3.3 Formation of Red Asymmetric Profiles for 
Molecular Spectral Line Transitions 

Why does our EECC dynamic model lead to the red profiles 
of molecular spectral lines? Before formal radiative trans- 
fer calculations, we would present a schematic analysis in 
Fig. [3] Line of Sight 1, which is close to the cloud core cen- 
tre, intercepts the ovals and straight arrow lines at points 
labelled 'B' and 'R' for the same approaching and receding 
projected velocity components, respectively. The intensity of 
an optically thick transition at a certain frequency will be 
dominated by emission from the nearest component to the 
observer, i.e. Ri and B3 along Line of Sight 1, because of 
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the absorption and scattering along the path. Assuming the 
local thermal equilibrium (LTE) for simplification, when the 
gas temperature increases towards the core centre (as for the 
case in the bottom panel of Fig. [2]), emission from point Ri 
will be stronger than that from point B3, which then gives 
rise to the asymmetric red profile feature in double-peaked 
molecular line profiles. On the other hand, as Line of Sight 2 
passes completely through the expanding envelope, no 'red 
profile' feature would be detected in this situation. 

For a more detailed analysis, we would first invoke the 
following radiative transfer equ ation (e.g. IChandrasekharl 
ll960l : lRvbicki fc LightmanHigTSH . namely 

= ~Iu + , (17) 

ar 

with r being the optical depth, S being the source function 
and I being the radiative intensity; subscript v denotes that 
these variables are all at a particular frequency v. Then the 
radiative intensity at blue and red shifted frequencies are 

/biue = 5i(l-e-^i)e"'"'+"'' (18) 
+5'2(l-e-"^)e-"='+5'3(l-e-"^) , 

/red = 53(1 -e""^)e"'"^+"^' (19) 
+S2(l-e-"^)e-"i +S'i(l-e-"i) , 

respectively, where subscripts 1, 2 and 3 denote variables 
at three respective positions along Line of Sight 1 shown in 
Fig. |3] Therefore the difference between the blue-shifted and 
red-shifted intensities /blue and /rod around this frequency is 

/blue - /red = (5*3 - Sl)[l + e" > ] (20) 

+ (S2-S'3)[e""^ -e-("i+"^'] 

Under the LTE condition and the Ray leigh- Jeans ap- 
proximation for a bla ckbody source function (e.g. 
iRvbicki fc Lightman|[l97^ . intensities can be expressed in 
units of brightness temperatures. Therefore the difference 
between the blue-shifted and red-shifted peak brightness 
temperatures Tbiue and T^cd is 

Tbiuc - Tred = (T3 - ri)[l + e-("i+"^+"3>] (21) 
+ (r2-T3)[e-"^ -e-("i+"^'] 
+ (Ti-r2)[e-"i -e-<"^+"^>] . 

We can immediately conclude from equation (21) that in 
an optically thick cloud environment (i.e. n S> 1, T2 S> 1, 
T-i ^ 1), the difference in peak brightness temperatures 
/blue — /red ~ T3 — Ti . Siuce the gas temperature gener- 
ally tends to increase towards the cloud centre, i.e. T2 > 
Ti > Ts, we then have Tbiuc — /rod < corresponding to 
the red-shifted peak intensity being greater than the blue- 
shifted peak intensity (i.e. an asymmetric red profile). In 
the opposite limit of optically thin regime (i.e. ri <g 1, 
T"2 ^ 1, T3 ^ 1), the difference in brightness temperatures 
/blue — /rod ~ 0, indicating that we will get equal peak inten- 
sities in blue- and red-shifted frequencies. When the blue- 
and red-shifted frequencies change (i.e. the loci of equal- 
projection velocity changes), all foregoing analytical results 
remain valid. Therefore we can apply these results to the 
entire blue- and red-shifted peaks in molecular spectral line 
profiles, and the following conclusion can be reached: double 
peaked molecular line profiles with stronger red peaks than 



blue peaks (i.e. red profiles) form under the optical thick 
condition in a dynamic EECC cloud core. If the outer ex- 
panding envelope is replaced by a static envelope, we would 
then have blue profiles for molecular transitions for optically 
thick lines along the Line of Sight 1 in Fig. [3] (see also Gao, 
Lou fc Wu 2009 for more details). 

4 DATA FITTING OF EECC SHOCK MODEL 

4.1 Molecular Spectral Line Profiles 

Adopting this EECC shock dynamic solution to model 
FeSt 1-457, we perform radiative transfer calculations for 
HCO+ (J= 1-0) line profiles u sing the RATRAN code 
jHogerheiide fc van der Takll20 00h with spherical symmetry. 
The Doppler b parameter in the code is taken as 0.1 km s~^, 
corresponding to an intrinsic line broadening of 0.17 km s~^. 
We introduced twelve shells with different thickness (more 
concentrated towards the cloud central region to account for 
the rapid variation of physical properties there) in perform- 
in g the calculations. The dust emissivity follows the model 
of lOssenkopf fc Hennin j (Il994h with bare ice mantles and 
with ~ 10^ years of growth. The abundance ratio of HCO"'' 
to H2 in this cloud core is listed in Table 1. Abundance vari- 
ations in radius r would affect molecular line profiles (e.g. 
iTsamis et a\\ \20oi\. Here a constant abundance with a cen- 
tral fractional drop (i.e. a st ep function) appears to be a sen- 
sible prescription (see also lEvans et al.ll2005l : iTafalla et all 
2006) . Physically, molecular abundance ratios vary with the 
initial metal composi tions and the evolu tion history of a 
molecular cloud (e.g. iTsamis et al.ll200^ : and the central 
abundance drop is caused by adhesions of molecules onto 
dust grjiiris in a relatively higher temperature environment 
(e.g. Tafalla et al.. 2006'). In the molecular spectral line pro- 
file fittings here, abundance values and the radii of abun- 
dance holes (see Table 1) in the step function serve as addi- 
tional fitting parameters to achieve reasonable agreements 
between our model results and the observational data. 

Given the beam width of ~ 24", we average our model 
spectral line profiles over a beam area of ~ 24" x 24" and 
multiply a beam efficiency of '-^ 0.75 for the IRAM 30m tele- 
scope (dotted curves in Fig. 1) to fit the line profile data. 
Fairly reasonable overall fits of model line profiles to those 
actually observed 'red profiles' (solid curves in Fig. [T]from 
lAguti et al.ll2007l ) support EECC shock dynamic structures 
(Fig. [2j for the dark globule FeSt 1-457. Red profiles and 
red-skewed single peak line profiles in Fig. [T] are encour- 
aging indicators of expansive motions present within this 
cloud core. In particular, the central 3x3 grids show good 
agreements between observations and the results of spheri- 
cal EECC shock model, suggesting that the spherical global 
expansion is at least a sensible first-order dynamic pattern 
characterizing the molecular cloud core FeSt 1-457. Certain 
deviations of spectral line profiles towards the right bot- 
tom grids in Fig. [1] then suggest that there might be a 
be amed outfiow or ce rtain rotations [as have been discussed 
in iGao fc Loul (|2010l )] overlapped with the global envelope 
expansive motions described by the spherically symmetric 
model. We may infer from the differences between model 
results and observations in Fig. [1] that this additional asym- 
metric motions are of the order of ~ 0.1 km s"'^. 



Three Diagnostics of EECC Star-forming Cloud Model 7 



Table 1. Molecular line transition properties: t'o is the rest- 
frame central frequency of the line transition; Xq is the molec- 
ular abundance with respect to H2 molecules in the cloud 
core FeSt 1-457; and rjj is the radius of the central abundance 
hole, within which a fractional abundance drop of 10^* is 
presumed. 



FeSt 1 -457 HCO"" (J=3-2) velocity 



Molec. 


Trans. 






Xo 




rh(cm) 


Note 


HCO+ 


J= 1 - 





89.19 


2x10" 


9 


7.5 xlO^" 




HCO+ 


J= 3- 


2 


267.56 


2x10" 


-9 


7.5 xlO^" 


b 


CS 


J= 2- 


1 


97.98 


3x10- 


-9 


1.05x10" 




CS 


J=3- 


2 


146.97 


3x10" 


-9 


1.05x10" 


b 




J= 1 - 





109.78 


1x10" 


-7 


6.0 xlO^^ 


c 



° Frequency valu es are taken from the LAMDA (e.g. 
ISchoier et al]|2005l') and abundance values are close to those 
of iTafaD^etah I l l2006h with slight variations. 
^ Theoretical model predictions for future observational tests. 

The abundance of C^*0 is enhanced to ten times the value 
of Xo for cloud radius r > 1.80 x 10^^ cm. 



FeSt 1-457 CS (J=2-1) 



O 
O 
m 



Velocity (km/s) 
) 5.5 6.0 6.5 



-24 
RA Offset (arc sec) 



-48 



Velocity (km/s) — 
5.0 5.5 6.0 6. 5 ^ 
4.0 




-24 
RA Offset (arc sec) 



Figure 4. Similar spatially-resolved grid-maps as shown in Fig. 
[1] but here for CS J= 2-1 (top panel) and Ci«0 J= 1-0 (bottom 
panel) molecular spectral line transitions. The receding velocity of 
cloud core FeSt 1-457 along the line of sight is estimated as 5.75 
km s~^. Spectral line fits are done using the same underlying 
polytropic EECC shock dynamic model described in the text, 
with step abundance profiles summarized in Table 1. 




24 -24 

RA offset (arc sec) 

FeSt 1-457 CS (J=3-2) 



Velocity (km/s) ^ 
5.0 5.5 6.0 6.5 ^ 



0.5 2 

Q. 

n 6 



RA offset (arc sec) 

Figure 5. Theoretical model predictions for the molecular transi- 
tions HCO+ J= 3 — 2 (top panel) and CS J = 3 — 2 (bottom panel) 
in spatially-resolved grid-maps of cloud core FeSt 1-457 on the ba- 
sis of our best-fit EECC hydrodynamic shock model (see Section 
3.2 for model details). Each grid has the velocity as the abscissa 
(with the line-of-sight receding velocity 5.75 km s"-"^ of FeSt 1-457) 
and the antenna temperature (for a 30m IRAM telescope beam 
efficiency of 0.75) as the ordinate. Model spectral profiles are con- 
volved with a typical beam width of ~ 24" . The spatially-resolved 
grid-map is centered at the right ascension aj2000 = 17''35™47.5'' 
and the declination i5j2000 = — 25°33'2.0", with a grid spacing of 
24". The predication of HCO+ J= 3 — 2 emissions cover a maxi- 
mum area of 168" X 168"; while for CS J = 3 — 2 emissions, the 
spatial distribution is 120" X 120". Molecular abundance infor- 
mation are summarized in Table 1. 



Meanwhile, spatially-resolved CS (J= 2-1) and C^^O 
(J= 1 — 0) spectral line profiles are also computed using 
the RATRAN code given the same underlying EECC shock 
dynamic model and fitted with observational data of FeSt 
1-457 shown in Fig. |4l The CS and C^**0 abundance ratios 
are also prescribed as step functions as listed in Table 1. In 
running the RATRAN code for radiative transfer calcula- 
tions, molecular transition data of HCO^, CS and C^*0 are 
all taken f rom the Leide n Atomic and Molecular Database 
(LAMDA; ISchoier et aU 2005). The reasonable fits of CS 
J= 2 — 1 and C^*0 J= 1-0 model molecular line profiles 
using the same EECC shock dynamic model to the observed 
corresponding line profiles (see Fig. |4)| confirm the viability 
of the underlying EECC shock model. In comparison with 
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spatially-resolved spectral line profiles of HCO"'" J= 1 — 0, 
the observed CS J= 2—1 and C^**0 J= 1-0 spectral line pro- 
files show less asymmetric features because of the relativel 
lower optical depths in these two molecular line transitions 
Both Figs. [1] and |4] show somewhat notable deviations be- 
tween model spectral line profiles and those observed at grids 
in the lower right of the map, which is most likely associated 
with non-spherical symmetry of any realistic astrophysical 
systems. We may infer from these complementary molecu- 
lar line model fittings, though with less asymmetries in line 
profiles, that the spherically symmetric EECC model repre- 
sents a reasonable description to this particular star-forming 
cloud. It is expected that more structural complexities on 
smaller scales (e.g. Redman et al. 2004 for rotations and 
Carolan et al. 2008 for bipolar outflows etc.) can complicate 
such radiative transfer fittings. 

As testable theoretical predictions for future telescope 
observations of the HCO"^ J= 3 — 2 and CS J= 3 — 2 emis- 
sion lines from the cloud core FeSt 1-457, we show our model 
calculation results as predictions in Fig. [5] Besides that all 
these predictions are based on the same EECC hydrody- 
namic shock model described in Section 3.2, abundance pro- 
files of HCO"^ and CS remain also the same as what we have 
used to calculate HCO"'' J= 1-0 and CS J= 2—1 molecular 
line transitions (see Table 1). Molecular transitions HCO^ 
J= 3 — 2 and CS J= 3 — 2 emissions have higher excitation 
temperatures than those for HCO^ J= 1-0 and CS J= 2—1 
emissions, respectively (e.g. ISchoier et al"]l2005l ). Therefore 
in the environment temperature of only ~ 10 K in a cloud, 
they are less excited and thus have lower emission intensities 
as shown in Fig. [5] The intensities of CS J= 3 — 2 emissions 
are even weaker because of a larger central abundance hole 
presumed (see Table 1). Because the abundance hole is large, 
less higher temperature cloud regions will contribute to the 
line emission, which makes the CS J= 3 — 2 emission very 
weak. This leads the conclusion that for certain molecules, 
the radius of the abundance hole will affect the intensity 
contrast between its transition lines with higher and lower 
excitation temperatures. We may also learn from Fig.[T]and 
Fig. [5] that optically thin molecular transition lines (here 
HCO^ J= 3 — 2) show less asymmetries in line profiles than 
optically thick transition lines (here HCO"'' J= 1 — 0). That 
is why optically thick lines are more sensitive diagnostics of 
hydrodynamic structures in molecular clouds. 



4.2 Column Density and Dust Continuum Profiles 

The dust extincti on observation of FeSt 1-457 of 
iKandori et al] l|2005h represents an important observational 
constraint to the model desc ription of this molecu lar cloud 
core. We infer from Fig. 2 of lKandori et all l|2005l ) that the 
cloud core FeSt 1-457 appears grossly spherical, yet with 
slight asymmetry towards the lower right corner. This asym- 
metry in the dust extinction observation coincides with the 
differences shown in spectral line profile model fittings as 
shown in Figs. [T] and U All these available observations 
might imply that there exists some kind of beamed outflows 



^ As discussed in iGao fc Loul 1 I2OIOI) . lower opacities would give 
rise to less asymmetries in molecular line profiles, thus not effec- 
tive indicators of cloud hydrodynamics in star-forming clouds. 



towards the lower right direction in FeSt 1-457, in addition 
to an almost spherically symmetric cloud core. As the de- 
viation of molecular lines in this part of the cloud core is 
toward red-shifted in the order of 0.1 km s~^ (see Figs. [T] 
and|4]), we estimate that a beamed slow outflow may recede 
from us with a typical speed of ^ 0.1 km s~^. 

As a further consistent model check, we compute the 
radial profile of column density using the same EECC shock 
dynamic model and compare it with that inferred from th e 
dust extinction data for FeSt 1-457 l|Kandori et all 120051 ). 
As the column density data has a spatial resolution of 30", 
our computed model profile is derived from the same EECC 
shock dynamic model with a convolution in 30" and shows 
a sensibly good fit (the top panel of Fig. [S| . As the radial 
profile of column density is sensitive to the outer radius i?out 
of the molecular cloud core and to the cloud number den- 
sity A*', this fit identifies a proper outer radius iiout and the 
molecular hydrogen H2 number density scale (see eqn (|16|) ). 
The outer radius of FeSt 1-457 used in the column density fit 
is i?out = 1.8 X 10'* an, which is also used in molecular spec- 
tral line profile fits and appears cons istent with that from 
the optical images (|Aguti et al.ll2007l ) at the same time. In 
short, the inferred variation trend of column density offers 
a constraint on a proper EECC shock dynamic phase with 
a grossly self-consistent radial density profile. 

In previous publications (e.g.lAlves. Lada fc Lada|[200ll : 
iHarvev et al.l l2003l : iKandori et al. 20051 ) regarding the col- 
umn density structure of molecular clouds, researchers nor- 
mally use a power-law radial density d istribution o r a static 
Bonnor-Ebert sphere (|Bonnoij Il956l : lEbertl Il955l ) etc. to 
describe density distributions within observed molecular 
clouds. These empirical/theoretical models also appear to 
provide reasonable fits to the observationally inferred ra- 
dial column density distributions, but none of them have 
been shown to successfully account for the observed molec- 
ular spectral line profiles from those molecular clouds, ex- 
cept that additional paramet erized conditions are intro- 
duced (e.g. iTafalla et al.ll2006t ). In contrast, our theoretical 
cloud models involving proper hydrodynamic flows and with 
self-consistent radial temperature va riations can give rise to 
asymmetric mole cular line profiles (|Gao. Lou fc Wul l2009l : 
ICao fc LorJ I2OI0I ) . According to the Bonnor-Ebert spher e 
model fits to a group of cloud cores in lKandorT et al.ll|2005l ). 
most of those best-fit static spheres are actually unstable 
Bonnor-Ebert spheres, implying that these molecular cloud 
cores are most likely undergoing hydrodynamic evolution in 
reality. Another data analysis bv lLi et"aLl l|2007l ) on massive 
quiescent cores in the Orion molecular cloud, gives a similar 
result that most of the cloud cores are unstable, if mod- 
elled as static Bonnor-Ebert spheres. All these results hint 
that we do need a theoretical model formulation for the hy- 
drodynamic evolution of cloud cores. That is, the spectral 
line profiles and dust emission data fitting are not incon- 
sistent with our EECC dynamic shock model. Essentially 
similar problems also exist in the model fits of dust con- 
tinuum emission observations described immediately below. 
Therefore to be consistently constrained by all available ob- 
servations, we explore the polytropic EECC hydrodynamic 
shock model here for the molecular cloud core FeSt 1-457. 

Radial profiles of (sub)millimetre continua at three 
wavelengths can be also computed using the RATRAN code 
to further check the consistency for the density and temper- 
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Figure 6. Radial profiles of column number density (top panel) and predicted (sub) millimetre continuum emissions (bottom panel) are 
shown here. The top panel compares the radial profile of column number density de rived from the EECC hydrodynamic shock model 
(solid curve) and the data with error bars inferred from dust extinction observations llKandori et al.ll2005l. dots with errors). Error bars 
are root-mean-square (RMS) scatters of data points within each annulus. The radial profile derived from the EECC hydrodynamic shock 
model is convolved with the observational spatial resolution of ~ 30". The same EECC hydrodynamic shock model that fits molecular 
spectral line profiles in Fig. 1 also offers a fair fitting (within RMS errors) to the radial profile of column density. The bottom panel shows 
continuum emission intensities at three wavelengths of 1.2 mm, 0.85 mm and 0.45 mm, which are predicted here for future observational 
tests. T hese (sub)millimetre continu um profiles are derived using the same underlying EECC shock dynamic model, with dust opacities 
given bv lOssenkopf fc Hennin^ ('l994) in the case of bare ice mantles and ~ 10^ yr growth, and convolved with a typical spatial resolution 
of ~ 10" X 10". Three pairs of gray curves along the three predicted continuum profiles show relevant ranges of uncertainties in the 
(sub) millimetre continuum emissions, which are estimated according to RMS errors in the column number density observation, using 
the relationship Iv oc p^'^. Abrupt turns around ~ 70" signal shock discontinuities and might be detected with high enough spatial 
resolutions. 



ature radial str uctures of clo ud cores. For these testable ob- 
servations (e.g. lShirlev et al.. 2000. 200^, for other cloud sys- 
tems) , we specifically predict such continuum radial profiles 
at three wavelengths of 1.2 mm, 0.85 mm and 0.45 mm us- 
ing the same underlying EECC hydrodynamic shock model 
for FeSt 1-457. Modelling such dust continuum emissions 
as optically thin with a constant dust opacity in r at a 
specified freque ncy v, the observed intensity at an impact 
parameter h is (|Adamslll99l1 : IShirlev et al1l2002l ) given by 

-2jy BAT,)j^^^£^^dr , (22) 

where Ta_{r) is the dust temperature taken to be the same as 
the H2 gas temperature T{r), p{r) is the mass density of H2 
molecules an d B,^(Td) is the Planck funct ion. The opacity 
K„ is that of lOssenkopf fc Hennin j (|l994l ) dust model for 
bare ice mantles with a growth time of ^ 10^ yr. As usual, a 
H2 gas to dust mass ratio 100 to 1 is adopted in our model 
computations. With the dust temperature and density pro- 



files from the EECC hydrodynamic shock model, we obtain 
radial profiles for 1.2 mm, 0.85 mm and 0.45mm continuum 
intensities. Results shown in the bottom panel of Fig. [S] are 
convolved with a telescope beam size of ~ 10" x 10" and 
a beam efficiency of 1. Error ranges of millimetre con- 
tinuum intensities are estimated by the root-mean-square 
(RMS) errors in the column density observation (the upper 
panel of Fig. |6]), through the relation oc p^'^. The abrupt 
'kink' at 70" from the centre reveals shock discontinu- 
ities of density and temperature in the EECC shock model 
(see the middle and bottom panels in Fig. [2]), and should be 
detectable with high enough spatial resolutions. 



5 CONCLUSION AND SPECULATIONS 

We have presented model results pertinent to three comple- 
mentary observational aspects towards the dark molecular 
globule FeSt 1-457; fairly reasonable fittings to asymmet- 
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ric 'red profiles' and to the column density radial profile 
derived from dust extinction data indicate that this dark 
globule most likely involve an early EECC shock phase of 
dynamic cloud evolution, while our predictions for spatially- 
resolved grid-maps for HCO^ J= 3 — 2 and CS J= 3 — 2 
respectively as well as the three predicted millimetre con- 
tinuum radial profiles can be specifically tested by future 
observations. The choice of our model parameter set is deter- 
mined by an extensive exploration and known quantitative 
as well as qualitative constraints. Slight adjustment around 
this chosen model parameter set is allowed for a refined im- 
provement of fitting. Although without a rigorous proof, it 
would be surprising that a dramatically different model pa- 
rameter set could reproduce all the available observational 
data so far. In addition to the spherical global EECC dy- 
namics, line profile and extinction observations might also 
suggest a slow beamed outflow receding from us at a typ- 
ical speed of ~ 0.1 km s~^ towards the lower right direc- 
tion in the grid-maps. The inferred underlying EECC hy- 
drodynamic shock model provides physical parameters for 
FeSt 1-457 cloud core. The central protostellar mass blob 
and the total mass of the cloud core are Mo ~ O.O2OM0 
and Mtot ~ 2.93M0, respectively; the central protostel- 
lar mass accretion rate is Mq ~ 0.69 x 10~^Mq yr~^ (see 
eq (O); the dynamical timescale of the cloud core FeSt 1- 
457 is estimated as t 2.8 x 10^ yr accordingly (see ex- 
pression (|16p V As an expanding shock travels through the 
molecular cloud core, abundance differences of certain tracer 
molecules (e.g. CH3O H and SiO) across the shock front (e.g. 
IJorgensen et al.ll2004l 'l could be detectable. 

As a strong supporting evidence to the present case of 
FeSt 1-457, a parallel theoretical model analysis on another 
star-forming molecular cloud cor e L1517B has also bee n con- 
ducted and completed recently (|Fu. Gao fc Loull2010l ). Sev- 
eral molecular line transitions of L1517B also manifest red 
asymmetric line profiles. On the basis of a shocked EECC 
similarity hydrodynamic model, several single point central 
molecular emission lines of HCO+(l - 0), HCO+(3 - 2), 
H"CO+(1-0), DCO+(3-2), H2CO(2i2-lii), H2CO(2ii- 
lio), CS(2 - 1), CS(3 - 2), SO(23 - 12), and SO(34 - 23) 
from the source L1517B [observed by 13.7 m telescope of 
the Five College Radio Astronomical Observatory (FCRAO) 
and TRAM 30 m telescope] are simultaneously fitted in 
an overall consistent manner. Taking into account of the 
fact that spatially-resolved HCO^ (J= 1 — 0) emission line 
profiles [observed by Delingha 13.7m telescope of the Pur- 
ple Mountain Observatory (PMO)], as well as 1.2 mm and 
850 ^m continuum emissions (observed by TRAM 30 m and 
JCMT respectively) can be consistently fitted by the same 
underlying EECC hydrodynamic shock model, we are fairly 
confident that the cloud core L1517B appears to be another 
cloud core that evolves in the EECC shock phase. Similar to 
what we have done for FeSt 1-457 here, predictions of sev- 
eral other spatially-resolved molecular transition lines, the 
radial profile of 450 /^m continuum emission, and visual ex- 
tinction (all based on the same EECC hydrodynamic shock 
model) are also presented for furthe r observational tests o f 
the model description of L1517B (see lFu. Gao fc Loul (|2010l ') 
for details). 

We speculate that other starless cores of grossly spher- 
ical molecular clouds with 'red profiles' might also involve 
EECC hydrodynamic shock phase in general. More sources 



need to be search and established by systematic model anal- 
ysis and data comparisons. The existence of global EECC 
hydrodynamic shock at early phases would affect the for- 
mation of protostars and their surroundings, including es- 
timates for star formation rates and initial mass functions 
that need to be further explored. 
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